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I. 

For a long time, it is generally believed that spin-spin interactions can only exist 
in a theory where Lorentz symmetry is gauged, and a theory with spin-spin in- 
teractions is not perturbatively renormalizable. But this is not true. By studying 
$h , the motion of a spinning particle in gravitational field, it is found that there exist 

spin-spin interactions in gauge theory of gravity. Its mechanism is that a spinning 
particle will generate gravitomagnetic field in space-time, and this gravitomagnetic 
field will interact with the spin of another particle, which will cause spin-spin in- 
teractions. So, spin-spin interactions are transmitted by gravitational field. The 
form of spin-spin interactions in post Newtonian approximations is deduced. This 
result can also be deduced from the Papapetrou equation. This kind of interaction 
will not affect the renormalizability of the theory. The spin-spin interactions will 
violate the weak equivalence principle, and the violation effects are detectable. An 
experiment is proposed to detect the effects of the violation of the weak equivalence 
principle. 
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1 Introduction 



It is know that the source of gravitational interactions is energy-momentum. In Newton's 
classical theory of gravity [1], the motion of a test particle in gravitational field is driven 
by the gravitational force on the particle, and the spin of the particle will not affect its 
motion in gravity. In general relativity [2, 3], the motion of a test particle in gravitational 
field is determined by geodesic equation. But the motion of a spinning particle is not given 
by geodesic equation. In general relativity, the motion of a spinning particle is given by 
the Papapetrou equation, which is deduced from the Bianchi identities and under the 
assumption of pole-dipole approximation [4, 5, 6, 7]. 

Quantum Gauge Theory of Gravity (QGTG) is first proposed in 2001 [8, 9, 10, 11]. 
It is a quantum theory of gravity proposed in the framework of quantum gauge field 
theory. In 2003, Quantum Gauge General Relativity(QGGR) is proposed in the frame- 
work of QGTG [12, 13, 14]. Unlike Einstein's general theory of relativity, the cornerstone 
of QGGR is the gauge principle, not the principle of equivalence, which will cause far- 
reaching influence to the theory of gravity. In QGGR, the field equation of gravitational 
gauge field is just the Einstein's field equation, so in classical level, we can set up its ge- 
ometrical formulation [15], and QGGR returns to Einstein's general relativity in classical 
level. The field equation of gravitational gauge field in QGGR is the same as Einstein's 
field equation in general relativity, so two equations have the same solutions, though math- 
ematical expressions of the two equations are completely different. For classical tests of 
gravity, QGGR gives out the same theoretical predictions as those of GR[16], and for 
non-relativistic problems, QGGR can return to Newton's classical theory of gravity [17]. 
Based on the coupling between the spin of a particle and gravitoelectromagnetic field, 
the equation of motion of spin can be obtained in QGGR. In post Newtonian approxi- 
mations, this equation of motion of spin gives out the same results as those of GR[18]. 
QGGR is a perturbatively renormalizable quantum theory, and based on it, quantum ef- 
fects of gravity[19, 20, 21, 22] and gravitational interactions of some basic quantum fields 
[23, 24] can be explored. Unification of fundamental interactions including gravity can 
be fulfilled in a semi-direct product gauge group[25, 26, 27, 28]. If we use the mass gen- 
eration mechanism which is proposed in literature [29, 30], we can propose a new theory 
on gravity which contains massive graviton and the introduction of massive graviton does 
not affect the strict local gravitational gauge symmetry of the action and does not affect 
the traditional long-range gravitational force[31]. The existence of massive graviton will 
help us to understand the possible origin of dark matter. 

The equivalence principle is generally believed to be one of the milestone of Einstein's 
general relativity. It has been tested by many experiments [32, 33, 34, 35, 36]. All these 
experiments give out null result for the violation of equivalence principle. On the other 
hand, it is generally believed that spin-spin interaction between rotating bodies will vio- 
late equivalence principle. Zhang et al. proposed a phenomenological model for spin-spin 
interactions [37]. Later, Zhang and his collaborator design a double free fall experiment 
to test equivalence principle with free falling gyroscopes [38, 39]. There has also been an 
interest in spin-gravitational coupling for a long time[40, 41, 42, 43, 21, 44]. 
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In a previous paper, the equation of motion of a spinning test particle in gravitational 
field is deduced based on the coupling between the spin of the particle and the gravito- 
magnetic field[44]. In this paper, we will study the post Newtonian approximation of that 
equation. The post Newtonian approximation of the Papapetrou equation is also studied. 
Then, we apply the equation to the problem of gyroscope which is moving around the 
earth. Finally, qualitative results on the motion of the gyroscope are given, and based 
on this calculation, an experiment is proposed to detect the effects of the violation of the 
weak equivalence principle. 



2 Equation of Motion of a Spinning Particle 

For the sake of integrity, we give a simple introduction to QGGR and introduce some 
basic notations which is used in this paper. Details on QGGR can be found in literatures 
[8, 9, 10, 11, 12, 13, 14, 16]. In gauge theory of gravity, the most fundamental quantity is 
gravitational gauge field C^(x), which is the gauge potential corresponding to gravitational 
gauge symmetry. Gauge field C^(x) is a vector in the corresponding Lie algebra, which 
is called gravitational Lie algebra. So C^(x) = C^(x)P a (fi, a = 0, 1, 2, 3), where C^(x) is 
the component field and P a = is the generator of global gravitational gauge group. 

The gravitational gauge covariant derivative is given by 

D li = d li -igC li (x)=G°d a , (2.1) 

where g is the gravitational coupling constant and matrix G = (G") = (5" — gC"). Its 
inverse matrix is G^ 1 = j} gC = {G~ ltl ). Using matrix G and G^ 1 , we can define two 
important composite operators 

g°f* = TrGZGZ, (2-2) 

gap = V^G-^Gp 1 ", (2.3) 

which are widely used in QGGR. In QGGR, space-time is always flat and space-time 
metric is always Minkowski metric, so g alS and g a p are no longer space-time metric. They 
are only two composite operators which consist of gravitational gauge field. The field 
strength of gravitational gauge field is defined by 

F» v (x) = ^-[D» , D v \=F^{x)-P a (2.4) 

where 

F; v = G%C: - G%C«. (2.5) 



In QGGR, gravitational gauge field is a spin-2 tensor field. The field equation of 
gravitational gauge field given by the least action principle is equivalent to the Einstein's 
field equation in General Relativity[12, 13, 14]. 
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According to literature [44], the equation of motion of a spinning test particle in 
gravitational field is 



Df 
~Dr~ 



= 9 



dx@ 
d7 : 



where J pa is the spin tensor of the particle, and 



Dr dr 



dr 



or equivalently 



Equation (2.6) can be written into the following form 

Dp' 



dx? 



Dt 



/a 
s i 



(2.6) 



(2.7) 



{21 



(2.9) 



where ff is the interaction force originated from the coupling of spin and gravitomagnetic 
field, 

ft = <T [fy (^sJ p °F s pa ) - d p {lg^J^F%^ dX " 



dr 



(2.10) 



3 Post Newtonian Approximation 

Now, let's discuss its the post Newtonian approximation. The standard procedure of the 
post Newtonian approximation can be found in text book[45, 46]. In gauge theory of 
gravity, post Newtonian approximation is done in a similar way. Let'e consider the case 
that a gyroscope moving around the earth. Then its post Newtonian approximations give 
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where and Q are given by the following relations 

= , 



C= -5" (a X J e ). 



(3.2) 
(3.3) 



In above relations, M e and J e are mass and spin angular momentum of the source. The 
leading terms of gravitational gauge fields are 



{ gCj = gq = -6 ij( /> + o(v*) 

gC° = + o(v 4 ) 

{ gCl = -gC° = \Q + o(v") 
The gravitational gauge field strength F" have the expansion 





5ik{dj<t>) - 


- 5 jk (di(j)) + oiv^/r) 




- \m 


- d,Q + o(v 5 /f) 


< 

gF? = 


-gF& = 


\diCk + S ik +o(v 5 /f) 




-9*Z = 


di<f> + o(v 4 /r) 



If we use the following relations, 



r>ot Tpa 



77a _ ™ 



we find that the gravitoelectromagnetic fields have the following approximations 

f = -\ £lkm (d m <P) + o(v A /f) 



gE» 



|(Vx () t + (v 5 /f) 
-(V0), + o(t) 4 /f) 



For spin tensor J' 1 ", we have 

5 



(3.4) 



(3.5) 



(3.6) 
(3.7) 



(3.8) 



(3.9) 



The quantity is defined by 



p7 1 7 5 ( ^9as_ %tf _ dg aP 
al3 2 9 \8xP dx a dx s J ' 

In post Newtonian approximations, its leading contributions are 

" Poo = §^ + d t (2^ + ^+C t +o(v 6 /r) 



1 jk 



1 00 



1 Qi 



r°. 



= \{djCi - did) - Sij +o(v 5 /r) 

= -Sijd k (f> - S ik dj(p + S jk di(p + o(v 4 /f) 

= (j) + i/j + C -V0 + o(v 7 /f) 

= drf + drt + oiv 6 /?) 

= -KdiQ + djQ-Sij't+o^lr) 



where ip is determined by the following equation 



(9 2 (A / 2 2 

2 - ^ +47iG[T 00 + T d 



(3.10) 



(3.11) 



<% 2 



(3.12) 



Using all above approximations, we can calculate the leading contribution of the force 
fs from equation (2.10), which is 



/; = d> 



^■(Vx c) 



J xV); 9 o 0+w-V0 +2j-(uxV)(9, 



-<9 O 0+ v -V0 



+ -(« -V) J-(Vx c) 



(3.13) 
(3.14) 



We can see that fl is of order v 3 J/r , and /° is of order v J/r . 



In general relativity, under pole-dipole approximation, the equation of motion of a 
spinning particle is given by the following Papapetrou equation 



^ (m^ + U/J — J - - J^R^ = 0, 



(3.15) 



where u a = ^ is the four-velocity, and -R^ 7 is the curvature tensor. Now, let's discuss 
its post Newtonian approximation. In most cases, 

DJ a P 



Dt 



< m, 



(3.16) 



so, 



Uf3- 



DJ al3 



<C mu a . 



(3.17) 



Therefore, in leading order approximation, the term up D £* 8 in the Papapetrou equation 
can be omitted. Then in post Newtonian approximation, equation (3.15) simplifies to 



Dp a 
~D7 



(3.18) 



Compare it with equation (2.9), we find that, in the present case, the interaction force 
originated from the coupling between the spin and gravitomagnetic field is 



After applying results (3.1) and (3.4), we can prove that 
'1 



J-(Vx c) 



=J .(„' X V) -d <f>+ v -V0 +^(u-V) J-(VxC) 



+ (J xV)i d (j)+ v -V0 +2 J .(u xV)(<9 4 0), 
1 



(3.19) 



(3.20) 



(3.21) 



Compare them with equations (3.13) and (3.14), we find that, in post Newtonian ap- 
proximation, the force f' s a is the same as ff. So, the equation of motion (2.6) and the 
Papapetrou equation give out the same results in the post Newtonian approximation. 



4 Spin-Spin Interactions Force 

Now, let's change the forms of equations (3.13) and (3.14) into more explicit forms. 
When we discuss the case that a gyroscope moves around the earth, the post Newtonian 

field and ( are given by equation (3.2) and (3.3) respectively. Then in leading term 
approximations, equation (3.13) is changed into 

SI = -3G(J- J ffi )^-3G(f ■ J ffi )i-3G(f ■ j)^ + 15G(f ■ J)(r • J e )^ 

+3GMe (l^k _ 3GM e (v ■ 7)^^ - 6GM ffi [j -(v x r)]£ + (j xV),f , 



and (3.14) is changed into 
/ s ° = -3G(v- J)(7- J ffi )4 -3G(y. J e )(r • jU- 3G(J- J e )(r ■ V)X 



(4.1) 



+ 15G(r • J)(r ■ J ffi )(f • v)^ - 3GM e [j -(v x r)](r • J •(« xV)f . 

(4.2) 
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In equation (4.1), the first four terms are from spin-spin interactions, the next two 
terms are from spin-orbit interactions. We can see that the spin-spin interaction force has 
the following form 

f„= —3G(J ■ J ffi )^ - 3G(r • J e )^ - 3G{r • j)^ + 15G(r • j)(r • J©)£, (4.3) 

which is symmetric under the exchange of two spin operators. So, the spin-spin interaction 
force is proportional to the spin of each particle, but decreases as the inverse quadruplicate 
of the distances. Equation (4.3) can be changed into another form 

S'„ = G^^L, (4.4) 

where 

M] k = -35 jk fi - 35ijf k - Zdikf-j + Ihfifjfk. (4.5) 

In above equation, f = is a unit vector. It can be seen that M l - k is symmetric under 
exchange of any two indexes of i, j and k. 



5 Test Gravity Theory With a Gyroscope 

Now, we discuss a classical test of gravity theory. Suppose that a gyroscope is placed in 
an orbit around the earth. For the earth, its gravitational field is quite stable, so the time 
derivative of the earth's gravitational potential almost vanishes 

In the test, we select a special orbit so that the spin axis of the earth is along the direction 
normal to the plane of the orbit, in this case 

r ■ J e =v ■ J ffi = 0. (5.2) 

In this case, the time component and the space component of the spin-dependent force 
fg and /] are respectively simplified into 



f s = -3G(J ■ JM-3G(7 ■ J)^ + 3GM ( 



(Jx«), 

) r 3 



(5.3) 



-3GM®(v ■ r)i^ii - 6GM e [j -(v x r0, 



fs = -3G( J ■ J e )(r • v)- - 3GM e [j -(v x 7)}(7 ■ v)\. (5.4) 

Then, we let the spin axis of the gyroscope be in the plane of the orbit. So, the spin axis 
of the earth J e is perpendicular to the spin axis of the gyroscope J, 

J® • J= 0, (5.5) 
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and r x v is perpendicular to the orbit plane, and therefore to the spin axis of the 
gyroscope 

j -(r x v) = 0. (5.6) 

In this case, the time component of the spin-dependent force /° vanishes, and space- 
component of the spin-dependent force finally becomes 

] = _3G(r • J)4^ + ZGmJ-I^P- - SGM^v • r)^p-. (5.7) 



where 



' A(r) 


= 1-20 


< £(r) 


= 1 + 20 


k C(r) 


2GJ® 
c 4 r 



The non-vanishing component of the inverse matrix of g a @ is 



9 tt = 
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B(r) 


grr = 


1 

A(r) 


g" = 


1 


A(r)r 2 


gfP = 


1 


A(r)r 2 sin 2 9 


g«* = 


<7* = ^ 



In post newtonian approximation, g a p is given by (3.1). In spherical coordinate system, 
we can prove that the non-vanishing components of g a p are 

' g tt = -B(r) 

g„ = A(r) 

gee = A(r)r 2 , (5.8) 

g w = A(r)r 2 sin 2 6 

, g<pt = gt<p = C(r) sin 2 9 

I I 2GM e 

r 

L _2GMe ( 59) 



(5.10) 



rTg is defined by (3.10). In spherical coordinate system, its non-vanishing components 
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are 
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(5.11) 



Now, let's discuss the equation of motion of the gyroscope, which is given by equation 
(2.9). Using above results on T^, we can obtain 



dV 

dr 2 



I A'(r) I dr\ 2 B'(r) /dt\ 2 
1_ 2 J 4(r) \At) "'" 2A(r) Ut / 



r + 



r 2 A'(r) 
2A(r) 



de\ 2 



(f) 



r + 



r 2 A'(r) 
2A(r) 



sin 2 9{^f -C'{r) sin 2 9 f r f r 



d 2 9 



C(r) n dt dip 
-4^ sin 20- — f- + 



dr 2 r 2 dr dr 

d 2 y? C'(r) dt dr 2C(r) 



'2 A'(r) 
- H — 

r A{r) 



drd9 . . n fd<py f e s 

- — : sin cos i-r-) = — , 

dr dr y dr J m 



(5.12) 



(5.13) 



dr 2 



+ 



r 2 dr dr 



dt d9 
dr dr 



A'(r) 
A{r) 



Ps 



dr dip „ d9 dtp 

— / + 2ctg£— / = ± 5.14 
dr dr dr dr m 



d 2 t B'(r) dt dr 3C(r) d 2 e dr d( P _ ft 
dr 2 B(r) dr dr r dr dr m 



(5.15) 



The spherical coordinate system is specially selected so that the polar angle 9 of the plane 
of the orbit is |. But the plane of the orbit can not always exactly be |, for there always 
exists small perturbation to the motion of the gyroscope. If the small perturbation to the 
polar angle of the plane of the orbit is denoted by SO, then we have 

9 = ^ + S9, (5.16) 

where S9 is a first order infinitesimal quantity. If we only keep first order infinitesimal 
quantity, then equations (5.12), (5.13), (5.14) and (5.15) are changed into 



d 2 r A'(r) 
+ 



dr 



dr 2 2A(r) \dr J 2A(r) \dr 



+ 



B'{r) (dt\ 2 {r 2 A{r))' fdip\ 2 C'{r)dp f ; 



2A(r) \dr J A(r) dr m' 



(5.17) 
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d 2 se | 2 dr dse | 5e /dfp\ 2 = f[ 

dr 2 r dr dr V dr / m ' 



dV 

dr 2 



+ 



"2 A'(r) 
_r A(r) 



dr dy d 
dr dr dr 



C(r) 



y/UM^LJ dr 



d 2 * 
dr 2 



+ 



£'(r) 3C(r)y/G'M @ L 
5(r) + ^ 



dt dr 
dr dr 



fl 
m 



(5.18) 
(5.19) 
(5.20) 



In order to go any further, we need to know the explicit form the the spin dependent 
force f s , which is given by equation (5.7). In the present coordinate system, the z axis is 
selected to be along the direction of J e , so 



J x v- 



eJ \l — - — smip sm{ip — ip ) + J cos(y? — (f ) 



r ■ J= rJcos(y? - ipo), 



(5.21) 
(5.22) 
(5.23) 



where tpo is the intersection angle between the spin J and semimajor axis, L is the semi- 
latus rectum, and e is the eccentricity. Then the spin-dependent force has the following 
form 

fl = 0, (5.24) 



where 



mi.}.} m;m .i^cm. i. 



c 2 r 4 



+ 



c 2 r 4 



cos(y? - <f )- 



(5.25) 
(5.26) 



The above expression for the force f s is in the Descartes coordinate system. We need to 

change its form into that in spherical coordinate system. The components of the force f s 
in spherical coordinate system are 



r fi = o, 
f e = -A 

J s r i 

{ ft = o. 

Then, equations (5.17), (5.18), (5.19) and (5.20) are changed into 



(5.27) 



dV A'{r) f 
dr 2 + 2A(r) \ 



dr\ 2 B'{r) SdtV (r 2 A(r))' fd(p\ 2 C'(r) dip_ 
d7y + 2A{r)\d^) 2A{r) \ dr" J ~ ~A{r) dr" ~ ' 



(5.28) 
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d 2 56 | 2drd69 | 5e = _f*_ 

dr 2 r dr dr V dr / mr ' 



dV 

dr 2 



+ 



2 A'(r) 
r A(r) 



dr dip d 



C(r) 



drdr dr y^GM^L J dr 



^ = 0, 



ft_ 

dr 2 



+ 



B'(r) 3C(r)y/0M 9 L 



B(r) 



+ 



dt dr 
dr dr 



0. 



(5.29) 
(5.30) 
(5.31) 



Equations (5.30), (5.31) and (5.28) yield three constants of motion, which are 

C(r) 



2*/ xdy? 
r Ayr)—— ■ exp 

dr 



„, ,dt 
B(r) — ■ exp 
dr 



y/GM @ L 



Jo (constant), 



c 4 r s 



1. 



(5.32) 
(5.33) 



Mr) ( ^V-4r+-r^-Te- 2C M/v / ^ 4G J @ J 2 _4O^0 = _^ (constant) 

(5.34) 



Multiply both sides of equation (5.29) with r 2 , we get 



dS9\ Jn 



dr V dr 



fa 



m 



Because 



d dip d J d 
dr dr dy? r 2 d</? ' 
we can change (5.35) into the following simpler form 



d<^ 2 



(5.35) 
(5.36) 

(5.37) 



The above equation (5.37) is the typical equation of motion of a forced oscillator. Applying 
(5.26), we get 

— | — = 6cos((/9 — <po)(l + ecosip), (5.38) 



with b a constant, which are defined by 



/ 3(CM e ) 3/2 3GJ e \ / J 

V C 2 J 2 LV2 C 2J2 L y ^ m 



(5.39) 



The general solution of oscillator equation (5.37) is 

59(ip) = cti cos ip + a 2 sin </? + «3(y?) cos + a^ip) sin 99, 



(5.40) 
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where a,\ and a 2 are two constants, and a 3 {p) and a^tp) are determined by the following 
two equations 

<ia 3 ((p) _ f s 



dp J^rri 
da 4 ((p) f s r 



simp, (5-41) 



cosy?. (5.42) 



dtp JqTti 
So, we have 

a 3 ((p) = b -^(p + ^^sm 2 (p-^^sm(pcos(p 

+ ebc^o (J _ cog 3 ^ + ebsm^, ^3 ^ 

a 4 (p) = - b -^p-*^ S m<pcos<p- b -^sm\ 

-e6cos^ (sin^- |sinV) - - cos 3 </?). 

Therefore, the general solution of equation (5.37) becomes 

59(p) = oil cos (p + a 2 sin p> — |y? sin(y? — <p ) — bsn ^° s in y> 
+y[-|cos</7 + cos(</7 + </?o) + |cos(2<^ - v?o)]- 



(5.43) 



(5.44) 



(5.45) 



The orbital period T of the gyroscope is 



T^/Jf-, (5.46) 
with clg the semimajor axis of the orbit of the gyroscope 

a G = ^j. (5.47) 

Suppose that at the starting point, (p — 0. Then at an arbitrary time t, the value of y? 
angle is about 

9irt 

pit) ~ — . (5.48) 

When the time i is many times larger than the orbital period T, we will have the following 
approximate solutions for a 3 and a 4 

a 3 (t) ~ ^ ^— fa sin v^o, (5.49 

a 4 (t) ~ - ^%^! ^cosy . (5.50) 
2a% 2 
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If at the starting point, we make the plane of the orbit exact at the plane with 9 — | , 
which requires that when t is small, or equivalently (p is small, 59 should almost vanish. 
Applying this initial condition to (5.40), we can obtain 



Ct\ = Ct2 = 0. 
Therefore, the solution (5.45) becomes 

66(t) = -9 m (t) sin(y?- ip ), 

with 

e m (t) 



3GM a 



3J e (GM a 



1/2 



2c 2 {La G yi 2 



2c 2 M S) L 2 aJ 2 



t. 



9 m (t) is the amplitude of the oscillation, which is linearly increased with time. 



(5.51) 



(5.52) 



(5.53) 



6 Observable Effects 



In the last chapter, we have discussed the motion of a gyroscope. When the plane the 
orbit is perpendicular to the spin axis of the earth, and the spin axis of the gyroscope lies 
in the orbit plane, the motion of the gyroscope in the eg direction obeys the forced oscil- 
lator equation of motion. The resonant force f s originates from the interaction between 
the spin of the gyroscope and the gravitomagnetic field of the earth. The solution of the 
equation is given by (5.52). Now, our question is that whether this forced oscillation is 
observable, or the amplitude 6 m (t) is large enough to be detectable? And which informa- 
tion can we obtained from this experiment? 

Inside the satellite that carries the gyroscope, the quantity 9 m (t) is not directly mea- 
surable, what is directly measurable inside the satellite is the relative displacement of the 
gyroscope inside the satellite. The relative displacement 5rg(t) is given by radius r times 
angular displacement 59(t), that is 



Therefore, we have 



with 



w m (t) 



r ■ 



Sr e (t) = r-59{t). 

Sr e (t) = -w m {t) sin(y? - tp Q ), 

3GM @ 3J ffi (GMe) 1/2 ' 
2c\La G fl 2 2c 2 M ffi L 2 a 3 J 2 



t. 



(6.1) 
(6.2) 

(6.3) 



In order that the oscillation can be observed, its amplitude w m should be large enough. 
Now, let's first estimate its magnitude. There are two terms in the right hand of (6.3), 
which represent two separate contributions. They are denoted by 



Wl (t) 



3GM e r 



2c 2 {La G y/ 2 \m 



J 



(6.4) 
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... 3rJa,(GM s )'/ 2 ( J\ 

WM = - 2^M 9 L^f UJ <6 - 5) 

where W\ represents the spin-orbit interactions, and w 2 represents the spin-spin interac- 
tions. So, we have 

w m (t)= Wl (t)+w 2 (t). (6.6) 
Suppose that the orbit is near the surface of the earth, in this case 

L ~ a G ~ r ~ ~ 6.38 x 10 6 m. (6.7) 

It is known that the spin angular momentum of the earth is about 

J e ~ 5.92 x 10 33 kg ■ m 2 /s. (6.8) 

Using all these approximations, we have the following estimations 

wAt) = 5.14 x 1(T 9 ( — J m/year, (6.9) 
\m J 

w 2 (t) = -1.01 x 10" 10 ( — j m/year, (6.10) 



where the unit of J/m is in m 2 /s. If J/m is about 10 3 m 2 /s, then 

Wi(t) = 5.14 fim/year, (6.11) 

w 2 (i) = -0.101 fim/year, (6.12) 
Both of them are large enough to be detectable. According to (6.6), we have 

w m (t) = 5.04 /im/year. (6.13) 

So, if the magnitude of J/m is about 10 3 m 2 /s, the oscillation amplitude of the gyroscope 
in the eg direction increases about 5.04 micrometer per year, which is detectable. 

In the above discussions, it is found that, the dominant contribution to the increasing 
of the oscillation amplitude comes from the spin-orbit contributions, which is about 51 
times larger than that of the spin-spin interactions. Now, the question is that, whether 
can we simultaneous detect the magnitude of both w± and u>2? The answer is that, if 
we want to detect both of them, we need at least two gyroscopes. Suppose that we have 
another gyroscope, which is in another satellite that is running around the earth in the 
opposite direction but with the same orbit parameters. In this case, the direction of the 
spin axis of the earth inverses, so the oscillation of the gyroscope in the ee direction is 

5r e (t) = -w' m (t)sm{<p-<po), (6.14) 

where 

w' m (t)=w 1 (t)-w 2 (t). (6.15) 
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After we detect both w m (t) and w' m (t), we can determine both spin-orbit contribution 
Wi(t) and spin-spin contribution W2{t) 

Wl (t) = W ^ +W '^\ ( 6 .1 6 ) 
w 2 (t) = Wm(t) : W - (t) . (6.17) 



7 Summary and Discussions 

In this paper, based on the equation of motion of a spinning particle in gravitational 
field, spin-dependent gravitational force is discussed. This kind of gravitational force con- 
tains two part, one is the spin-orbit interactions, and another is spin-spin interactions. 
Both of them are spin-dependent force. This force originates from the coupling between 
the spin of the particle and gravitoelectromagnetic field. So, whether the results in this 
paper are reliable depends on whether the coupling between spin and gravitoelectromag- 
netic field is correct. It is known that, another important effects of gravity, the preces- 
sion of a gyroscope in gravitational field, also originate from the spin-gravitomagnetic 
coupling[45, 46, 18]. That is, the influence of the spin-gravitomagnetic coupling to the 
motion of the spin is the precession of a gyroscope, while its influence to the motion of 
particle is what we discussed in this paper. It is know that the effect on the precession of 
a gyroscope had already been detected by experiments recently [47, 48], and experimental 
results are quite consistent with theoretical expectations, which tells us that the spin- 
gravitomagnetic coupling does exist and is reliable. As long as the spin-gravitomagnetic 
coupling has influences to the motion of the spin, it must have influence to the motion of 
the particle. 

The spin-dependent force is expressed by equations (3.13) and (3.14). The existence 
of this force explicitly showed that not all gravitational actions are inertial actions, which 
means the violation of weak equivalence principle [49]. We can consider this violation from 
a traditional way. We may consider that the satellite which carries the gyroscope is a free 
falling frame of reference. According to the weak equivalence principle, no objects in 
the free falling frame of reference can fell gravitational force, and therefore all objects in 
that frame should be relatively static or uniform rectilinear moving. But now, we know 
that, owing to the spin dependent force, the gyroscope is neither relatively static nor 
uniform rectilinear moving, but oscillating with increasing amplitude. The existence of 
this oscillation means that the gyroscope feels gravitational force in the free falling frame 
of reference, which violates the weak equivalence principle. 

For a long time, it is generally believed that, even if there exist spin-spin interaction, 
it must be too weak to be detectable. Indeed, the spin-spin or spin-orbit interaction force 
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is extremely weak. According to equation (5.7), the spin-spin interaction force is about 



3GJJ ff 



If we suppose that the angular momentum of the gyroscope is about lO^kgm 2 / s, then 
this force is about 8 x 10~ 17 iV, which is indeed too weak to be directly detectable, and is 
many orders smaller than the upper limit given by experiments [32, 33, 34, 35, 36]. But, 
this force will cause a resonance oscillation in the eg direction, which is detectable. Direct 
measurement of this oscillation will be of great importance, for it is the first direct obser- 
vation of the violation of the weak equivalence principle. There is a long standing dispute 
on the reliability of the weak equivalence principle. Some treated it as a corner stone of 
general relativity, while some others criticized that it is not strictly hold[49, 50, 51]. Now, 
the experiment can serve as a direct experimental solution of the dispute. 



Because QGGR can return to general relativity in classical level, and the corner stone 
of QGGR is gauge principle, our conclusion is that the equivalence principle is not the 
necessary starting point of general relativity. Therefore, the violation of weak equivalence 
principle does not affect the correctness of general relativity. On the contrary, if equiva- 
lence principle is still considered as corner stone of general relativity, the theory will be 
logically self-contradictory, for in that case, the violation of weak equivalence principle is 
a result of the equivalence principle. (Because Papapetrou equation is a direct result of 
Einstein's field equation.) 
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